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tency with covariant gauge Feynman rules. 

11.15Bt,12.38Bx 



1 



Typeset using REVT^X 



I. INTRODUCTION 



Perturbation theory for non-Abelian gauge fields is usually formulated in covariant 
gauges. There the Feynman rules nearly can be read off from the momentum space rep- 
resentation of the Lagrange density including gauge fixing and Faddeev- Popov ghost terms 
(c.f. any textbook on gauge field theory, e.g. [Q). The only subtlety is the proper pole 
prescription for the propagator 1/k 2 which is achieved by the causal e prescription, i.e. by 
the physically motivated requirement that particles (antiparticles) propagate forward (back- 
ward) in time. 

Despite the great success of covariant gauge Feynman rules there are situations where 
a different gauge choice is preferable. However it turned out during the last 15 years that 
things are much more involved if perturbation theory ought to be formulated in non-covariant 
gauges. The most prominent example is the case of temporal gauge where an additional pole 
(1//cq) appears in the gauge field propagator. Assuming the absence of ghost contributions 
and temporal propagator modes it was shown first in [[J that the actual value of gauge 
invariant quantities depends on the pole prescription for this additional propagator pole 
(see also || and references therein). According to j| this pole is caused by the fact that 
the temporal gauge condition is insufficient to completely fix the gauge degrees of freedom. 
A complete gauge fixing condition, on the other hand, automatically regulates this "gauge 
pole", albeit breaks translational f| or rotational || invariance, which at least makes the 
use of temporal gauge Feynman rules less attractive. 

A second example is the Coulomb gauge where singularities appear in the energy in- 
tegrations of two loop diagrams @-@|. These singularities must be cured by a very subtle 
prescription |8||| or alternatively by the introduction of new multi-gluon vertices [[U]]. As 
already pointed out in [jTOf these problems are due to operator ordering ambiguities which 
are difficult to handle in the path integral approach. 

These examples show that it is far from being trivial to find consistent sets of Feyn- 
man rules for non-covariant gauges. With "consistent" we simply mean that a given set of 
Feynman rules reproduces the same value for any gauge invariant quantity as the covari- 
ant gauge Feynman rules. The consistency of different sets of momentum space Feynman 
rules was systematically studied by Cheng and Tsai BJ7|, pT|1 . A diagrammatic scheme was 
developed in P,|TT| which serves to prove without analytic calculations that different sets 
of Feynman rules yield the same results for gauge invariant quantities. Cheng and Tsai 
called this the "theorem of equivalence" . To apply this scheme to a given set of Feynman 
rules the propagators must not have any singularity, i.e. every singularity must be regulated 
and the regulator must be kept finite during the course of calculation. E.g. for the case of 
temporal gauge it was shown |7| that some prescriptions of the gauge pole mentioned above 
reintroduce ghost fields and/or temporal propagator modes which vanish if the regulator is 
taken to zero but nevertheless contribute to loop integrals as long as the regulator is kept 
finite. Hence the discrepancy between covariant and temporal gauge calculations observed 
in H is due to the neglect of ghost and temporal propagator modes in the temporal gauge 
using a specific gauge pole prescription, namely the principal value prescription. 

In this article we will present a set of Feynman rules for the radial (Fock-Schwinger) 
gauge and apply the diagrammatic scheme of Cheng and Tsai to this set of rules to show its 
consistency with covariant gauge Feynman rules. The radial gauge with its gauge condition 
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X 



»A"(x) = 



(1.1) 



was introduced long time ago [0,|l3|] and rediscovered several times (e.g. Hl4j). It found 



widespread use in the context of QCD sum-rules (e.g. []I5| ) since it enormously simplifies 
the task of organizing the operator product expansion of QCD n-point functions in terms 
of gauge invariant quantities by expressing the gauge potential via the gauge covariant field 
strength tensor. 

In ( |1 . 1|) the origin is singled out by the gauge condition. Therefore this gauge condition 
breaks translational invariance and perturbation theory cannot be formulated in momentum 
space as usual but must be set up in coordinate space. For this reason only few efforts have 
been made to establish perturbation theory for radial gauge |16|-pl|] . Most of the discussion 



was concentrated on the question about the properties of the free radial gauge propagator 
while less efforts where made to explore the other Feynman rules, especially to figure out 
if ghost fields truly decouple. Recently a radial propagator with surprising properties was 
proposed in pi]] . 

0. 



In Sec. Ill we 



In the next Section we briefly recapitulate the results found in 
generalize the diagrammatic scheme of Cheng and Tsai, originally developed for momentum 
space, to coordinate space Feynman rules. In Sec. [Iy| we apply this scheme to a complete 
set of radial gauge Feynman rules including the radial gauge field propagator proposed in 
An additional check for the consistency of this set is presented in Sec. [V| using the 
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Slavnov- Taylor identities of radial gauge. Finally we summarize our results in Sec. |V|. 



II. THE RADIAL GAUGE PROPAGATOR 

In the following we work in a D-dimensional Minkowski space, i.e. with one time and 
D — 1 space dimensions. 



In |2T|] it was shown that the (full) radial gauge field propagator can be expressed in 



terms of a gauge invariant Wilson loop in the following way: 

K^V) ■■= (A;(x)A b u (y)) = S a \l N _ 1( \ 2 Jim dfdlW^x'^y') . (2.1) 

where TV is the numbers of colors and the Wilson loop W\ connects the sequence of points 
0, y', 0, x', x, along straight lines (c.f. Fig. [j]). As a Wilson loop with cusps and self- 
intersections W\ has so-called "cusp singularities" (c.f. [^2] and references therein) appearing 
at quadratic (and higher) order in the coupling constant. Due to the prefactor 1/g 2 in ( ^T|) 
these singularities show up at order g°, i.e. already for the free radial propagator. The 
explicit form of the latter is given by 



(A;(x)Al(y)) 

g^u[{x -yf -ie] 



T(D/2-l) b( 2 _ ,_ D/2 



4tt d / 2 
i 



d; Jdsx u \{sx - y) 2 - ie] 1 - ' 2 - &t Jdty, [(x - tyf - ie} 1 - 1 " 2 
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dffi Jds Jdtx >y[(sx- tyf -ie] 1 ' ' 2 ^ . (2.2) 
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A-D 

which is obviously singular for D —> 4. That the free radial propagator diverges in four 
dimensional space was already pointed out earlier in W7\. The reason for this singularity 



was presented in [21]. It is a surprising result that the renormalization properties of a special 



class of Wilson loops influence even the free radial propagator. However it was argued in 



21[| that this singularity is even necessary to reproduce correct results for gauge invariant 
quantities calculated in radial gauge. 

Note that we have used dimensional regularization for convenience. Of course the cusp 
singularities of Wilson loops and therefore also the singularity of the radial propagator might 
be regularized in a different way as well. In any case the regulator has to be kept finite during 
the calculation which however is required for loop calculations anyway. In addition a finite 
regulator is necessary to make the diagrammatic scheme of Cheng and Tsai applicable to 
radial gauge Feynman rules. 



III. THE THEOREM OF EQUIVALENCE IN COORDINATE SPACE 

The theorem of equivalence was originally formulated by Cheng and Tsai for momentum 
space Feynman rules 0. In this section we generalize it to coordinate space. The basic 
statement is that all gauge invariant quantities are independent of the arbitrary function 
A M (x, y) if the following Feynman rules are used: 
The free gauge field propagator 

y) = sab K>> y) - ^ A ^ y) - d l A ^ x )] > t 3 - 1 ) 

the ghost field 

G? c (y, x) = -igf ahc [{g, v U x - dffiWiy, x) + d v y D^(x, y)} , (3.2) 

and the common three and four gluon vertices. Here denotes the Feynman propagator 
given in Appendix the "ghost field" is the product of ghost-gluon vertex and ghost 
propagator, and the coordinate space representation of the three and four gluon vertices is 
also given in Appendix [A|. Note that we get the expressions for Feynman gauge for vanishing 
A M . Hence the ghost field ( |372|) can be written as 



Gf c (y, x) = -tgf abc (g, u n x - fl^)A"(y, x) + G* abc (y, x) . (3.3) 

To prove that all gauge invariant quantities are independent of A M we observe that in the 
expression for the gauge propagator (|3.1| ) A M is always accompanied by a partial derivative. 
This reflects the fact that the gauge modes and the longitudinal modes coincide in the 
free field case of a gauge theory. By partial integration we arrive at an expression where 
this partial derivative acts on a vertex which is connected to the propagator. The partial 
derivative acts on a three gluon vertex (Kl) as 
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dt T$l(w, x, y) = igf abc [5(w - y)(g^n w - d™d;)5{w - x) 

- S(w - x)(g^n w - d™d:)S(w - y)] . (3.4) 

Obviously the operators inside the brackets project on transverse components, i.e. 

(g^n w -d:d™)d» = 0. (3.5) 

These operators act on the next propagator which is connected to the vertex at x, /x, b or 
y ) c, respectively, e.g. 

igf abc jd D x (g Ufl B w - d:d;)5{w - x)D$ (x, x') 

= -gf ab ' c g/K™ - J) + dx> Gf c (x f , w) . (3.6) 

This yields ghost type structures as well as contractions —gf a y c gy8{w — x'). As we shall 
show in the following the former cancel gauge dependent contributions from ghost loops 
while the latter are related to four gluon vertices. From ( |3^4|) and flcTBD we get 



^ d D xd D y d>T&(w,x,y)D%\x,x')D™' (y,y>) = 
Jd D y [-gf ab ' c g/S(w - x') S(w -y) + dg Gf c (x>, w) S(w - y)} D%'(y, y') 
- Jd D x [- g r bc 'g/S(w - y') 5(w -x) + d»: G? c '(y\ w) S(w - x)} D$ (x, x') . (3.7) 

So far this looks very complicated. However a graphical representation of this formula 
is rather simple and indeed all calculations which are necessary to prove the theorem of 
equivalence can be performed graphically. Once all rules are stated, there is no need to 
write down any analytical expression. Formula ( |3.7|) is depicted in Fig. In the following 
all lines labeled with D denote gluon propagators. Three and four gluon vertices are dotted. 
G denotes a ghost field which includes a vertex and a propagator line according to the 
definition ( fT2| ). An arrow at the end of a propagator line denotes a partial derivative acting 
on the respective vertex. Finally an arrow appearing as a label denotes the contraction 
mentioned above, i.e. the first contribution on the r.h.s. of ( p.6|) . 

Instead of presenting a complete set of diagrammatic rules first and applying it to an 
example afterwards we will choose a simple example right now and introduce the necessary 
rules step by step as soon as we need them. As an example we take the two loop vacuum 
contributions shown in Fig. |3|. The combinatorial factors for these diagrams are 1/12 for the 
diagram involving three gluon vertices, 1/8 for the four gluon vertex diagram and 1/2 for 
the diagram with a ghost loop. It is useful to state the relative statistical weights and signs 
of all diagrams explicitly. This is depicted in Fig. [4|. An overall factor 1/12 is suppressed. 
In addition the diagrams are symmetrized with respect to the propagator lines. 

To prove the gauge invariance of the sum shown in Fig. |4| we have to show that the result 
remains the same if all D's and G's are replaced by D F and G F , respectively. To do that we 
will even prove a somewhat more rigorous statement: The sum remains the same if all fields 
D and G labeled by the same number are changed to D F and G F no matter of the gauge of 
the fields with different numbers. Thus we will prove in the following that without changing 
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the actual value of the whole sum D\ and G\ can be replaced by Df and Gf without any 
reference to the explicit form of D 23 and G 2 ^. This is enough to prove gauge invariance 
since this scheme can be applied to the diagrams of Fig. |] consecutively for three times: In 
the first step D\ and G\ are changed to Df and Gf. Next D 2 and G 2 can be replaced by 
D 2 and G 2 in the same way now dealing with Df , Gf , D3, and G 3 . Finally L> 3 and G 3 are 
replaced by their Feynman gauge expressions starting from diagrams where the lines labeled 
by 1 and 2 already are in Feynman gauge. Thus we will prove that the gauge of one line in 
all diagrams can be changed without changing the value of the sum of diagrams. 

We start with the first diagram of Fig. || and decompose D\ in its Feynman and A parts 
according to ( |3.1|) . This is shown in Fig. [5] where a partial integration has already been 
performed changing the relative sign of the A contributions with respect to the Feynman 
part. 

Now we use the diagrammatic rule of Fig. || to obtain Fig. |B|. 

In this way we get two types of diagrams. The first type involves ghosts and partial 
derivatives (arrows). The latter act on the second three gluon vertex and we can use ( |3.4|) 
again. In addition to (|3.6| ) we need a relation where the propagator D is replaced by A. 
However this is nothing but ( |3.3|) . Using (|3^ ), (|3.6| ) and (|3^3|) we get 



JdPx d D y d x w Tfc(w, s, y)A"(a/, x)D${y, y f ) = 

Jd D y [G^ abc (x\ w) 5(w -y)- G?*{x>, w) 6{w - y)} D%(y, y') 

- jd D x [-gf abc 'g/S(w - y r ) S(w -x) + d»y, Gf* (y\ w) S(w - x)] A»(x\ x) (3.8) 

which is shown in Fig. 0. 

For a diagram of Fig. § with a ghost line we get Fig. As a first example how the 
cancelation of gauge dependent diagrams shows up we note that the second diagram on the 
r.h.s. of Fig. || is identical to diagram G of Fig. |. These diagrams cancel each other (note 
the minus sign in front of the third diagram of Fig. |6|). 

The last diagram of Fig. |8| can be transformed in a diagram with a ghost loop in the 
following way: The action of the arrow (partial derivative) on the ghost vertex can be split 
up in a sum of two terms using the product rule of differentiation.^ On the one hand the 
partial derivative acts on a ghost field as (cf. flOl) ) 



d^G^iy, x) = d£G» Fabe (y, x) = -gf abc S(x - y) , (3.9) 

and on the other hand it acts on the upper line. This is shown in Fig. || The corresponding 
equation is 

Jd D xS(v - x) [G^ bc (y, x)5(x- w)] = 

- Jd D xS(x-w)d;[G^ c (y,x)S(v-x)}- Jd D xS(v - x) gf abc S(x - y) S(x - w) (3.10) 



x It is interesting to mention that the momentum space equivalent of this product rule is the 
energy-momentum conservation at the vertex. 
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and the symbol 1 in Fig. |9| denotes the last integrand. 

In this way the last diagram of Fig. |8| yields two new diagrams as shown in Fig. [U]. Since 
the only difference between 1 and the contraction arrow is a g^-term it is easy to get the 
last line of Fig. [U]. 

Now we collect all diagrams which can be deduced from the two diagrams involving G"s 
shown in Fig. [6| and compare the result, Fig. [II, with the original diagrams of Fig. [|: In 



addition to the ten diagrams of Fig. [TT| we get ten analogous diagrams not shown here from 
the last diagram of Fig. |. We refer to these diagrams as the "reflected" diagrams of the 
ones shown in Fig. [IT]. Diagrams a and b of Fig. [TT] replace G\ by Gf in Fig. [|G. The same 
is achieved by f and g in Fig. |I. The corresponding reflected diagrams transform Fig. f|F 
and J in the same way. Diagrams c and j as well as e and h cancel each other. Diagram 
d and the reflected diagram of i replace D\ by Df in Fig. |4|H. In the same way Fig. |]E is 
transformed by Fig. |ll|i and the reflected diagram of Fig. [TTld. 

Thus for the diagrams E- J of Fig. [4] we have already succeeded in replacing D\ and G\ 
by D[ and Gf . We are left with the four gluon vertex diagrams B and C on the one hand| 
and with the ghost-free diagrams of Fig. ^| and their reflected diagrams on the other hand. 
In Fig. |12| we have decomposed D\ of Fig. [4|B using ( |3.1|) . 

Obviously we need a relation which connects contraction arrows with four gluon vertices. 
Indeed the following relation holds 

- jdPw' gf abe g\S(v - w) S(v - w r ) T^{w\ x, y) 

+ JdPx' gfaecg^>5(v - X) 5(V - x') T^\x\ V , w) 

+ Jd D y f gfaedg\>S(v - y) S(v - y) T^(y\ w, x) = d° K Q K J^{v, w, x, y) . (3.11) 

This is shown in Fig. [13|. Connecting the upper legs with each other and also the lower legs 
we find that the last but one diagram of Fig. ^| and its reflected diagram exactly cancel the 
last diagram of Fig. [12].(] The same holds true for the diagrams where D 2 is replaced by D^ >) 
i.e. Fig. [4p instead of B and the second diagram of Fig. [6| instead of the fourth. 

We conclude that the last two diagrams of Fig. |5| can be decomposed in a way that the 
resulting diagrams can be used to replace D\ and G\ by Df and Gf in Fig. |}B-J. This 
proves the theorem of equivalence for our example of two loop vacuum bubbles. 

In the same way the theorem of equivalence can be proven for higher loop vacuum 
diagrams, scattering amplitudes and Wilson loops. There are a few more diagrammatic rules 
collected in Appendix [B] which we did not need here for our explicit example. However they 
become important for higher order Feynman diagrams. The complete set of diagrammatic 
rules, Figs. [2], |7|, |^, [13|, [15], [16|, and [H], can be applied to any gauge invariant quantity in any 



order of the coupling constant. The necessary steps to prove the gauge invariance of a sum 
of loop diagrams are as follows: 



2 Diagram D remains unchanged since no D\ appears there. 



3 Note that the third diagram of Fig. 13 vanishes due to its color structure if the upper legs are 
connected. 
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Symmetrize all diagrams with respect to all lines (as shown in Fig. [|). 

Start with the diagrams which have only three-gluon vertices. Decompose one gluon 
propagator according to ( [Tip . 



• Use Fig. ^ and [7| as long as there are <9 M -arrows pointing at a three-gluon vertex. 

• Use Fig. y to pull arrows out of closed ghost loops (c.f. Fig. |10|). 

• At this stage all A M -dependent diagrams which have neither contraction arrows nor 
four-point vertices should drop out. 

• Finally, diagrams with different topology are related by the rules of Figs. |13|, [15], [16|. 

To get an idea which diagrams are related it is instructive to contract the lines with 
contraction arrows. 



For Wilson loops the additional relation shown in Fig. [17] has to be taken into account. 



This scheme can also be generalized to cases where the Yang-Mills fields are coupled to 
matter fields. In principle for each new vertex two new diagrammatic rules must be found, 
which have to show how the two introduced arrows (denoting the partial derivative and 
the contraction) act on this new vertex. For the case of full QCD (albeit using momentum 
space Feynman rules) this complete set of diagrammatic rules is given in [[23|] and applied to 
three loop diagrams. Independently of the original work of Cheng and Tsai || a similar set 
of diagrammatic rules was developed recently by Feng and Lam also for momentum space 
Feynman rules and applied to QCD |24| and electro- weak theory |25| to simplify perturbative 



calculations by choosing an appropriate gauge. 

With the diagrammatic scheme developed in this Section we return to the case of radial 
gauge Feynman rules. 



IV. APPLICATION TO RADIAL GAUGE 

Now we have developed all tools necessary to present a consistent set of radial gauge 
Feynman rules. The free radial gauge field propagator is given in ( |2.2|) . To make the scheme 
of Cheng and Tsai applicable the regulator (here D — 4) has to be kept finite until the end 
of the calculation. Note however that the regularized propagator already fulfills the radial 
gauge condition. This is in contrast to the case of temporal gauge. There the temporal mode 
of the propagator vanishes only if the regulator vanishes. This mode, however, couples to 
the ghost fields. As we shall see in the following ghost fields decouple for radial gauge since 
already the regularized propagator strictly satisfies the gauge condition. 

Comparing ( p.l|) with (|2.2| ) we get 



i 11 
A u (x, y) = x a Jds D F av (sx, y) - -d« Jds Jdtx-y D*„(sx, ty) . (4.1) 



After some simple algebraic manipulations we find that according to ( p^l) the appropriate 
ghost field is given by 
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GT(x,y) = -igf 
It is easy to see that 



x v Jds8(sx-y)-d^D^(0,y) 



(4.2) 



G„ c i x , v) ~ Vv (4.3) 

holds. On account of the gauge condition ( |1 . 1|) we find 

G? c (x,y)DZ'(y,y') = 0. (4.4) 

Thus the ghost fields decouple from the gauge fields. Of course one expects that to happen 
for gauge fields obeying an algebraic gauge condition. As mentioned in the introduction 
this is however a non-trivial problem for the case of temporal gauge. This example teaches 
us that it is not sufficient to simply study the Faddeev- Popov ghost determinant. A more 
careful examination, e.g. like the one presented here, is necessary to state the decoupling of 
ghost fields. 

One subtlety has to be clarified before we can apply the scheme developed in the last 
Section to our case at hand: One frequently used tool for the manipulations of Feynman 
diagrams was the partial integration of derivatives d*. For all the statements of the last 
Section to be valid we have to make sure that this partial integration produces no surface 
terms. Clearly this depends on the actual form of A. Contrary to other radial gauge 
propagators presented earlier the propagator (|2^ ) is well-defined at the origin (for a detailed 
discussion of this point c.f. |2lH). Thus the behavior at the origin cannot be a source of 
surface terms. However things are more complicated at infinity. A careful inspection of 
the propagator (|2^ ) shows that for large x the first and third contribution decreases like 
l/(x 2 ) D / 2-1 . Also the propagator for the Feynman gauge (|A2| ) shows this rate of decrease and 
there no problems with surface terms appear. However the second and fourth contribution 
of (|2.2| ) drops only like l/(x 2 Y D ~ 3 ^ 2 . This weaker rate of decrease is always accompanied 
by a partial derivative d*. The worst case that might appear is that all legs at a vertex are 
of this type. A typical example for such a coincidence is 

'd D w jd D x Jd D yT^(w,x,y)d^Ay(w,w')d^(x,x')dyA ul (y,y'). (4.5) 



Inserting (|^|) and evaluating the 5-functions we are left with a D dimensional integral and 
an integrand with a partial derivative and three d A contributions. Hence the behavior for 
large x is given by Jdx/(x 2 ) D ~ 7 ^ 2 which might cause a logarithmic surface term for D — > 4. 
Before discussing this case in more detail we note that this actually is the only problematic 
case. If any of the OA contributions is replaced by a term which falls off like l/(x 2 ) D / 2_1 no 
surface term shows up if D is in the vicinity of 4. If there is any surface term it would be 
caused by a partial integration in (|45| ). Fortunately we can calculate this potential surface 
term explicitly by inserting the definition for the three gluon vertex (|A4|) and using ( |3.4| ). 
We find 

fd D w jd D x Jd D yT^(w,x,y)d™A x ,(w,w^ 

+ JdPw jd D x (4.6) 
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Thus also in this problematic case no surface term appears. 

We conclude that for the manipulations of radial gauge Feynman diagrams partial in- 
tegrations can be performed without any problems. Therefore the diagrammatic scheme 
presented in the last Section can be applied to the radial gauge Feynman rules. This proves 
that a consistent set of Feynman rules for the radial gauge is given by the radial propagator 
( ^2|) and the usual three and four gluon vertices and (\KEj), respectively. The ghost 

fields decouple (for the temporal gauge this depends explicitly on the pole prescription) 
and no additional vertices are necessary (for the Coulomb gauge one has to encounter that 
problem) . 

V. SLAVNOV-TAYLOR IDENTITIES 

Slavnov- Taylor identities for Yang-Mills fields formulated in radial gauge have already 
been derived in [JT6| . Besides the diagrammatic scheme presented in Sec. |T| these identities 



also serve to check the consistency of a set of Feynman rules. In addition these identities are 
an important ingredient to set up a renormalization scheme for radial gauge perturbation 
theory. For this latter aspect we refer to the remarks made in [|16[ . 



In the following we only study the consequences of the Slavnov- Taylor identities for the 
gluon self energy. In principle also higher one-particle irreducible n-point functions can be 
studied. We start with the generating functional for Yang-Mills fields in radial gauge and 
introduce a Lagrange multiplier field C: 

W[J] := jvA 5{x^{x)) 

x exp (i jd D x (£ym[A](x) + j;(x)A^x))^j 

= Jv[A, C] exp (i jd D x (C YM [A](x) + C a (x)x,A»(x) + J£(x)A£(x))) . (5.1) 

In addition we introduce the generating functional for vanishing source terms, 

W :=W[J = 0], (5.2) 
and full and free Green's functions for an arbitrary quantity £: 

<£) := ^ fv[A,C\Zexp{iS[A,C\) (5.3) 

and 

(0o := (0 

respectively. Here the action S is defined as 

S[A,C] := ld D x(C YM [A](x)+C a (x)x fl A^x)). (5.5) 
It is well-known that the operator 



(5.4) 

9=0 
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0^(x) = g^U x - d^8 u x (5.6) 

which appears in the free part of £qcd is n °t invertible due to the gauge invariance of the 
theory. However in the enlarged space of A and C fields the operator O can be inverted. 
We denote propagator and self energy in the larger space with 



(C a (x) Al{y)) (C a (x)C b (y)) 



_ D*{x,y) (A;(x)C b (y))\ 
{(C%x)Al(y)) {C%x)C b (y))) 

and fl@ $ (x,?/), respectively. The large Greek letters denote numbers from to 4. For to 
3 they coincide with the Lorentz indices of the gauge fields while 4 denotes the additional 
Lagrange multiplier field C. 

The Dyson-Schwinger equation turns out to be 

D = D + D UD , (5.8) 

where "multiplication" includes summation over all indices and integration over the space 
time variables, i.e. 115 is an abbreviation for 

4 N 2 -l 

E E d D xTl<g s (w,x)Df c *(x, y ). (5.9) 

E=0 6=1 

D is the free contribution of D. The inverse is the operator which appears in the free action 
(cf. (ED): 

(ZV)^ (x, y) = -i5 ab ( °^ X) x * ) S(x - y) . (5.10) 



The Dyson-Schwinger equation fl5.8|) can be transformed to 

n = Do 1 (D - Do)D- 1 , (5.11) 
which can be further evaluated using the Slavnov- Taylor identities 

{A«(x)C b (y)) = (A;(x)C b (y)) (5.12) 

and 

{C a (x)C b (y)) = (5.13) 

stated in ]16[. We find 

D-D =( D - D °° ); (5.14) 



hence: 
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n = -i 




(O(D-D ) 0\ =_! _ / 0(D - Do) 
1 { x(D - D ) ) ~ 1 { 



D~ 




(5.15) 

where we have used the gauge condition ( |LT1) in the second line. Stars denote arbitrary 
values. Since the self energy is symmetric with respect to an exchange of all variables we 
conclude that it has a non-vanishing contribution in the gauge sector, only: 

n=fnoy (5 , 6) 



In addition we get the transversality of the self energy by multiplying ( |5.15| ) with (<9 M ,0) 
from the left 

(011,0) = (<9,0)H) ( 0{D ~ Do) ° j D- 1 = (0,0) , (5.17) 



due todO = (cf. (p|)). 

We conclude that as a consequence of the Slavnov- Taylor identities the radial gauge 
gluon self energy obeys the transversality condition 

d*Il^(x,y) =0 (5.18) 

We will finish this Section by checking this equation for the lowest order contribution in the 



coupling constant. Fortunately we can use the diagrammatic rules of Sec. |TJ. Since the 
ghost fields decouple only two diagrams contribute to the self energy at order g 2 . They are 
shown in Fig. [14]. To derive the r.h.s. we have made use of the fact that the first diagram 
of Fig. [F| is symmetric with respect to an exchange of upper and lower line. In addition 
we have applied the diagrammatic rule of Fig. 0. Using the relation shown in Fig. O we 



find that the last diagram in Fig. [F| is exactly canceled by the last but one diagram. The 



remaining diagram involving G vanishes since G is contracted with a radial propagator D, 



i.e. 



jd D y Gf c (x, y) Dg{y, y') - y v D£(y, y')=0. (5.19) 



This proves that the g 2 contribution to the gluon self energy indeed satisfies fl5.18[) . Again 



this demonstrates the consistency of our set of radial gauge Feynman rules as well as the 



usefulness of the diagrammatic rules presented in Sec. Ill 



VI. CONCLUSIONS 

In this article we have presented a complete set of radial gauge Feynman rules for non- 
Abelian gauge fields. The rather surprising results concerning the radial gauge field propa- 
gator have already been discussed in ETJ : 
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• The propagator can be expressed in terms of a gauge invariant Wilson loop; 



• it is singular in four dimensional space; 

• this singularity is necessary to reproduce correct results for gauge invariant quantities. 
The other properties of this set of rules are: 

• Feynman rules must be formulated in coordinate space since the gauge condition ex- 
plicitly breaks translational invariance; 

• ghost fields decouple; 

• besides the conventional three and four gluon vertices no additional vertices are nec- 
essary. 

Introducing such a new set of rules it is crucial to show its consistency with the conventional 
covariant gauge Feynman rules, i.e. to prove that results for gauge invariant quantities 
calculated with both sets of rules agree with each other. To achieve that we have generalized 
the theorem of equivalence to coordinate space Feynman rules. Originally it was developed 
by Cheng and Tsai to check the consistency of different sets of Feynman rules formulated 
in momentum space. 

To the best of our knowledge such a systematic study of coordinate space Feynman rules 
has never been presented before. In |20| the consistency of radial gauge Feynman rules was 



checked by an explicit one loop calculation using basically the same technique as presented 
here. However it was not pointed out there how this technique could be generalized to other 
Feynman diagrams involving higher order loops. In addition the radial gauge propagator 
used in was later shown to be incorrect |T7| . 



A few additional remarks concerning the singularity of the radial gauge propagator ( ^2|) 
are necessary: In principle there are three sources for a possible singular behaviour. First, the 



already mentioned cusp singularity discussed in great detail in pT|| . Second, an additional 
singularity shows up if x and y are aligned with respect to the origin. These divergences are 
already present in Euclidean space. In Minkowski space a third type of singularity appears 
if x — y is light-like. Of course it would be interesting to study also the latter two types of 
singularities in more detail. At least the third one is also connected to Wilson loops with part 



of the contour on the light cone [^] . Despite the lack of a more detailed examination of these 



additional singularities we note that all these singularities are regularized by dimensional 
regularization. Thus also these divergences do not spoil the arguments presented in this 
article. 

To set up perturbation theory in radial gauge a complete set of Feynman rules is of 
course only the first step. Next one must figure out how renormalization works in this 
gauge. Due to the lack of translational invariance one might suspect that local counter 
terms are insufficient to renormalize the theory. However some arguments that this is not 
the case are given in . Further general considerations as well as explicit loop calculations 



of self energies and vertex corrections are necessary to work out a renormalization scheme 
for radial gauge. We expect that the diagrammatic scheme presented here might serve to 
simplify such calculations. 
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APPENDIX A: FEYNMAN RULES IN COORDINATE SPACE 



Since Feynman rules are commonly used in momentum space we have collected some 
coordinate space expressions for propagators and vertices in this Appendix. 

In Minkowski space with one time and D — 1 space dimensions the propagator of a scalar 
particle is given by p7 



d D k 



-ik-(x—y) 



(2tt) d k 2 + irj 



Y{D/2- 1) 



[{x - yf 



16 



l-D/2 



47F D/2 

The Feynman propagator got additional factors of unity in Lorentz and color space 

(B^x)Bl(y)) Feyn = S ab D^(x, y) = S ab g, v A s (x - y) . 
We define the (full) three gluon vertex T as 
(A\{w) A%x) AZ{y)) = 

fd D w' d D x> d D y> (A a x (w) AtWHAfr) A^WHWv) ^V)>tSvK v') I 



(Al) 



(A2) 



(A3) 



the four gluon vertex Q is defined in a completely analogous way. 

The expressions in lowest order of the coupling constant (bare vertices) are given by 



T^(w, x, y) = -igf ahc [g x ^ w - d v x ) 5(x - y) 5(w - y) 

+ g^(d* - c£) S(y - w) S(x - w) 
+ g ^-d^5(w-x)8(y-x) 



(A4) 



and 



Qjt%(v, w, x, y) = -tg 2 \f abe fc d e(9 Kfl 9 uX ~ g KV g X ») 



+fa d efUg K V ~ 9^9" X )} 5(v - w) 5(v - x) S(v - y) . (A5) 
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APPENDIX B: ADDITIONAL GRAPHICAL RULES 



Here we complete the set of diagrammatic rules presented in Sec. |T|. In addition to 
Fig. ^, [7], ^, and [T3| we have to discuss the case that a contraction arrow points at a four 
gluon vertex (Fig. [15]) or at a ghost gluon vertex (Fig. p~6f) . The corresponding formulae are 



d D v' gfeah9 K K ^(u -v)6(u- v') Qhb X c7( V '> w , x ,y) 
+ Jd D w' gf ebh g x x ,S(u - w) 5{u - w') Q K ^{v, w \ x > v) 

= Jd D x' gfehcg^n'Siu -x)5{u- x') Q^ v (v, w, x', y) 
+ jd D y' gfehdg u v>5(u - y) 5(u - y') Q K ^h\v, w, x, y') 



(Bl) 



for Fig. 15 and 



d D v' Gf e (w, v') 5(v - v') gf ecd S(x - v') S(x - y) = 
d D v' Gt be (w, v') 5(y - v') gf eca g\8{x - v) 5(x - v') 
+ jd D v> G? c (w, x) 5(v> - x) gf eda g\8(y - v) S(y - v') . 



(B2) 



for Fig. [II]. The last relation is a somewhat complicated way to write down the Bianchi 
identity for the structure constants.^ It is missing in 0. 

For Wilson loops we have to encounter the case that a partial derivative (arrow) "meets" 
a contour integral, e.g. 

ill 

Jdcfx Jda 2 Jda 3 Q(a 1 > a 2 > cr 3 ) w^ 2 (a 2 ) w^ 3 (a 3 ) 



x Jd D Xl Jd D x 2 Jd D x 3 D^(w(a 1 ),x 1 )^ a2) A^(w{a 2 ),x 2 )D^(w{a 3 ),x 3 ^ 
x T%%£{x 1 ,x 2 , x 3 ) Sp(t ai t a2 t a3 ) 

1 cri cri 

= Jdcix Jda 3 Jda 2 w^{ai)w^ 2 {a 2 )w^{a 3 ) 

cr 3 

X Jd D Xl Jd D X 2 l^^^^^^l)^!)^)^ 2 ^^),^)^ 3 ^^),^) 

x T^( Xl ,x 2 , x 3 ) $v(t ai t a2 t a3 ) . (B3) 

Here w(a), a £ [0, 1], parameterizes an arbitrary Wilson loop contour. We make use of the 
relation 



(B4) 



4 jabe jcde _j_ jace jdbe _|_ jade j 



ace -pdbe 



cade -pbce 



o. 
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to evaluate the contour integral: 



1 u\ 

Jda x JdazW^G^w^Gz) Jd D x x Jd D x 2 Jd D x 3 D fIlU1 (w(a 1 ),x 1 ) 



o o 



x [A^(w(a 1 ),x 2 ) - A^(w(a 3 ),x 2 )} D^(w(a 3 ),x 3 ) 

X ^vlvlv^ ( x l ) X 2i x 3) Sp(tai^02^a 3 ) • 



This is depicted in Fig. |TJ. In addition it is shown how the result is related to other diag 
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FIGURES 

FIG. 1. The Wilson loop W\ which is connected to the radial gauge field propagator. 
FIG. 2. Partial derivative acting on a three gluon vertex. 
FIG. 3. Two loop vacuum diagrams. 
FIG. 4. The sum of vacuum diagrams which ought to be gauge invariant. 
FIG. 5. Graphical decomposition of a propagator. 
FIG. 6. Application of the diagrammatic rule shown in Fig. |2|. 
FIG. 7. How the partial derivative acts on a three gluon vertex with A M on one leg. 
FIG. 8. Application of the graphical relation shown in Fig. ^. 
FIG. 9. Product rule of differentiation at a ghost gluon vertex. 
FIG. 10. How to pull an arrow out of a closed ghost loop. 
FIG. 11. Calculation of all G-contributions from Fig. [| A. 
FIG. 12. Decomposition of a propagator in a four gluon vertex diagram. 
FIG. 13. A relation connecting three and four gluon vertices. 
FIG. 14. Transversality of the one loop self energy contribution. 
FIG. 15. A relation connecting four gluon vertices. 
FIG. 16. A relation connecting ghost vertices. 
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FIG. 17. Additional relation for Wilson loops. 
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